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ABSTRACT 

We study how the effects of quantum corrections lead to notions of irreversibil- 
ity and clustering in quantum field theory. In particular, we consider the virtual 
"charge" distribution generated by quantum corrections and adopt for it a statistical 
interpretation. Then, this virtual charge is shown to (a) describe a system where the 
equilibrium state is at its classical limit {h ^ 0), (b) give rise to spatial diffusion of the 
virtual cloud that decays as the classical limit is approached and (c) lead to a scenario 
where clustering takes place due to quantum dynamics, and a natural transition from 
a "fractal" to a homogeneous regime occurs as distances increase. 
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Mercader and W. Zurek, Cambridge University Press (1993). 
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In quantum field theory, the dynamics stems from both the interactions among the 
fields and the quantum fluctuations to which they are subject. Quantum fluctuations 
are perhaps the most fundamental feature of quantum field theory, affecting fields, 
their sources and the vacuum in which they evolve. 

Induced by the fact that quantum fluctuations are to a certain extent random , 
we will study how this affects irreversibility. Namely, since true randomness impairs 
one's ability to carry out an exact reconstruction of the past history of the system, 
it is legitimate to expect some irreversibility at the microscopic (quantum field) level 
because of the presence of quantum fluctuations in this domain. 

One possible way to study this irreversibility is to introduce some fundamental 
statistical notions into the realm of the quantum field theory. For example, one can 
consider the effect that quantum fluctuations have on the "charge" density^ of the 
Poisson equation satisfied by the quantum corrected potentials. The quantum fluctu- 
ations affect the interaction energy, and this induces a modified "charge" density; this 
is how the Uehling potential appears in QED. The statistical notions come in when 
we interpret the induced charge density as a probability density, an action that we are 
justified in taking due to the mathematical properties of the charge density induced 
by quantum-fluctuations. We associate with the full interaction energy a "charge" 
density which (as is done in cosmology and astrophysicsB) we interpret as a prob- 
ability density characterizing the spatial distribution of virtual charges surrounding 
the original charges. By doing this, we are making use of the rich conceptual depth 
of quantum field theory and explicitly taking into account that, in the computation 
of the quantum corrections at some scale, we are actually including not only simple 
two body processes with a fixed and precise impact parameter of the order of that 
scale, but the effects of many body interactions with impact parameter smaller than 
the distance scale that we are probing. Because of this, divergences appear in the 
theory which renormalization removes and re-interprets. The physical potential that 

one measures is thus subject to fluctuations, and it is impossible to specify the exact 

•^We enclose the word charge in quotes because we refer to it in a generic sense, that is, we are 
not specificaUy referring to electric charge. We have in mind the charge for the source of the Poisson 
equation satisfied by the effective interaction energy. 
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dependence of the potential on the individual components of the virtual "charge" 
density. We may, however, ask about the probability of occurrence of some particular 
configuration or similar questions of a statistical nature. This goal is achieved by 
adopting for the virtual "charge" charge density a statistical interpretation in term of 
a probability distribution whose density is the quantum corrected "charge" density. 
By studying the effective (in the RG sense) form of the interaction energy, one 
acquires information on how the properties of the virtual cloud change with scale 
and, through the statistical interpretation, gain insight on how changes of scale affect 
the dynamics of the cloud in aspects relating to the approach to equilibrium or even 
structure and form generation properties. 

We can obtain the effective interaction energy by solving its RGE. The solution to 
this equationS yields the interaction energy as a function of the scale at which one 
probes the system, and contains the modifications due to the presence of quantum 
fluctuations. Unfortunately, the corrections that are included to all orders of pertur- 
bation theory are only the leading logarithms; but, as is well known, even with this 
limitation, many interesting physical consequences can still be extracted. Since en- 
ergy has canonical dimension of inverse time and no anomalous dimension, the RGE 
has the solution. 



Here V is the interaction energy, tq is the distance between the interacting sources, 
and a is a reference distance. The quantity ^o(-^) is the effective coupling and A the 
scale parameter; go{X) satisfies the RGE Xdgo/dX = —P{go)- 

In general, for a massless mediating field (such as photons or gluons), the effective 
interaction energy for two point particles separated by a distance a is given by 



Computing to 1-loop order, where /3 = Podh choosing X = r/a, and taking a to 



Strictly speaking, ct is a function of the distance at which the system is probed. This can be 
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V{r,go,a) = C^a~'^r 



-l+a 



(3) 



a is related to the /5-function for the couphng go, and to one-loop is given hj a = 
— 2/5o5'o- three essential properties: (i) it is proportional to h, (ii) it vanishes 

smoothly in the limit of zero coupling and {Hi) as defined, it is a function of the 
momenta which essentially counts the number of degrees of freedom excited from the 
vacuum at distances less than tq. 

For the case of QED one recognizes here the first term in the short distance 
expansion of the famous Uehling potential; in QCD one recognizes the expression for 
the interquark potential. 

According to potential theoryB, Equation @) satisfies a Poisson equation whose 
right hand side is proportional to the "charge" density dictated by quantum correc- 
tions. This charge density describes how the phenomenon of vacuum polarization has 
modified the vacuum and given rise to the formation of virtual pairs that, as is well 
known, affect the strength and properties of the interaction in a substantial way. 

However, we have no information as to the actual distribution of these pairs, but 
by interpreting the vacuum polarization charge density as a "probability density", 
we can derive information on the virtual cloud and its physics as a many body (sta- 
tistical) system. From the mathematical point of view this is possible thanks to the 
intimate relationship© that exists between potential theory and probability theory. 
The relationship between the potential and the charge density (away from r = 0) is 
through the Poisson equation = -|-47rp(r), which for our isotropic potential gives 



where A' is a constant. We only need to require from the density p{r) that it be 
a positive and integrable function on its support in order that it can actually be 
interpreted as a probability density (see below). 

made explicit by computing it in a mass dependent substraction procedure'!! . However, it will be 
sufficient for our present purposes to take it to be a constant. 




(4) 
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We now examine some properties of Eq.(^. First of all we see that p(r) is the 
solution to the functional equation 



p(Ar) = A^p(r) 



(5) 



with j3 = —3 + cr; this implies right away that p(r) describes a fractal distribution 
of charge embedded in 3-dimensional configuration space, and with a Hausdorff (or 
fractal) dimension given hj df = +a . This is not surprising, since a has its origin in 
the deviation from canonical scaling due to the quantum fluctuations. Furthermore, 
since a changes as we change the size of the domain on which we probe the virtual 
cloud, it turns out then that the Hausdorff dimension also changes and we are then 
dealing with a multifractal. For QCD and Quantum Gravity in its asymptotically 
free regime, a is positive. In QED and other non-asymptotically free theories, a is 
negative. When the Hausdorff dimension is positive, one sees intuitively that there is a 
natural tendency to suppress these fluctuations since, contrary to what happens in the 
case of a negative Hausdorff dimension, the configuration space "cannot accommodate 
them" : it is not big enough! 

To reveal some of the physical consequences of the randomness associated with 
these quantum fluctuations, we will study a few of the features of the "statistical me- 
chanics" of the density in Eq.(^). In particular, we will look at the entropy associated 
with p and also will give a glance at the nature of the stochastic processes that it 
supports. 

Requiring normalizability of the probability density, we get 



p = Ar 



-3+0- 



(6) 



with 



A 



47r 



'0 



(7) 



for cr > 0, and 
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when a < 0. Here Rq denotes an IR cutoff, necessary in the case of positive a, in order 
to ensure the finiteness of the probabihty distribution. For non-asymptotically free 
theories, vq is an UV cutoff which becomes necessary for the same reasons. The IR 
cutoff can be identified, e. g., with a typical hadronic size and ro with Planck's length. 
The resulting probability density is of the Pareto typeE^; this is a natural consequence 
of the renormalization group origin for p(r), which is ultimately responsible for the 
functional equation in and the associated scaling. Scaling leads to Levy-type 
distributions, which turn Pareto in some limit. Notice also that both densities go to 
zero in the classical {h 0) limit. 

To gain information about the "equilibrium configurations" supported by these 
distributions, we construct and compute an entropy-like quantity. We will assume 
that, as in thermodynamics, the state of equilibrium corresponds to the maximum 
for the entropy. In other words, we will assume that it reveals a preferred "direction" 
for stability in the evolution of the quantum field system. It is possible to discuss two 
types of entropies: a "coarse grained" entropy based on the probability distribution, 
and a "fine grained" entropy or "differential" entropyQ, based on the probability 
density p(r). We will discuss only the latter, which is defined through 



with the integral extended over the full support of the variable r. The constant C 
needs to be introduced because p(r) is a dimensional quantity, and we do not have 
the equivalent of a Nernst theorem to set a reference valid for all physical systems. 
This constant, as in thermodynamics, fixes the minimum entropy of the system. 
We will set Boltzmann's constant k equal to one. 

The differential entropy for the case when a is positive gives 
^In the same sense as in information theory^. 




(9) 
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and for negative a one gets, 

5("<°) = 1 log^C + 31ogro (11) 
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From these two expressions we see that the "entropy constant" , C, may be related 
in an interesting and useful way to the "extreme" volume (the volume at the cutoff) of 
the quantum system. Such a choice has the advantage of eliminating the dependence 
of the entropy on the cutoff. For asymptotically free theories, if we choose C pro- 
portional to the maximum volume, i.e., proportional to the volume of the IR-cutoff 
to the cube, then the dependence on the cutoff and extreme system size disappears 
from the expression for the entropy. The equivalent statement is also true for the 
case of non-asymptotically free theories, where the "Nernst-volume" corresponds to 
the minimum volume that can be physically reachable. 

These entropies are shown in Figures 1 and 2. 

We see that in asymptotically free theories, the equilibrium state corresponds to 
the largest possible value|^ of a. Now, since matter fluctuations contribute negatively 
to 0", the most stable state corresponds to a configuration (or system size) where the 
matter states do not contribute. Because of the decoupling theoremlI9, this occurs 
for the maximum-possible size of the system: that is for its classical limit! 

On the other hand, when the theory is non-asymptotically free, the maximum of 
the entropy happens when a goes to from the left. This means either the classical 
limit, in the sense that h 0, or that one probes the system at distances so large 
that no quantum fluctuation contributes to a. 

Thus, for both, asymptotically-free and non-asymptotically-free theories one sees 
that the maximum of the entropy is attained at the classical limit. 

How does this transition to the equilibrium state take place? Since diffusion is a 
mechanism for the transition to the equilibrium state for a system where randomness 

^For positive a, S^"^^^^ has a maximum at cr = 3; unfortunately, this value lies outside the range 
where we can reasonably trust the approximations made here. 
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is at work, we will examine diffusion in quantum field theory, and check if the emerging 
picture is consistent with what we have obtained from the entropy. In quantum field 
theory, the fractal nature of the probability density leads us to expect some form 
of fractal diffusion or, more precisely, diffusion through fractal brownian motion. 
Because of this "fractality" we expect the diffusion to be non-space filling, and to 
lead to either clustering {a > 0) or screening when the system relaxes; it is clear that 
this behavior ought to depend on the sign of a. In what follows we will examine some 
of the properties of the coefficient of diffusion derived from p(r); then we will examine 
the correlation integral in some specific cases. 

The virtual charge density cloud described by the probability density Equation 



(|^) will "diffuse" in a manner similarllll to what occurs in Brownian motion. We 
expect that the components of the virtual cloud scatter randomly among themselves; 
the collision probability density is p(r); the virtual charges execute random walks 
controlled by p(r); the net effect of these processes is that the virtual cloud diffuses 
outside of the ball of radius R. 



As is well knownllli the probability that a "particle" executing a random walk, 
with each individual step governed by a probabihty density p(r), be after A^-steps at 
a position between r and r + df, is given by 

W{r)£f= [AnDtr^^^ exp [- \f\'^ / (ADt)] £f (12) 
where D is the coefficient of diffusion, 

n is the number of steps per unit time and N = nt is the total number of steps, (f^) 



is the second moment of the probability density p(r). The function W{r) in Eq. (|1^) 
satisfies the 3-dimensional diffusion equation 



dW n 



dt 
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with boundary condition that iy| = at an absorbing medium, or vanishing normal 



derivative of at a reflecting surface: ^ 



0. How far and how fast a disturbance 



propagates depends on the coefficient of diffusion: for large D one has large amounts 
of diffusion, and viceversa. When = 0, there is no diffusion. 

In addition to the the Markovian nature of the process, there are two assumptions 



involved in the derivation of (JT^) which we must bear in mind: (a) the number of steps 
N is very large and (6) in the computation of the characteristic function, it is assumed 
that kr « 1. These two assumptions mean in our case that the time intervals that 
we consider are "long" compared with virtual time intervals and that the size of the 
region in which we expect the diffusion to take place, is "large" compared to the 
Compton wavelength about which we are computing the coefficient a. Both of these 
restrictions are amply met by the validity of our approximations. 

In order to compute the coefficient of diffusion on needs to obtain the second 
moment of the probability density. From a radial distance Sq to a radial distance 
s > So, the second moment of p(r) is 



2 + a 

and therefore, in some cases, the coefficient of diffusion diverges (naively) as s goes 
to infinity. This again is not surprising, since our probability distribution is the limit 
of a Levy-type distribution, and thus its moments are divergent. 

As previously, we classify the situation depending on whether we have an asymp- 
totically free theory or not. When a > 0, we can take sq = and set s = Ro, the IR 
cutoff. The coefficient of difussion is given by 

^(<x>0) ^ !^^Rl (14) 

In non-asymptotically free theories, we take sq = tq (the UV-cutofF) and leave s 
unspecified. Then, 



62 + (T 



S \ " (Tq^"^ 



(15) 
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The expression for Z}('^>'^) reflects the divergence in the second moment of Levy- 
type distributions. However, quantum corrections again play an important role here: 
physics limits the size of the system {Rq cannot go to infinity) and in the classical 
limit, when a —>■ 0, difussion stops. Before this happens, i.e., when the typical energies 
are larger than the inverse Compton wavelength of the IR-cutoff , the system will tend 
to diffuse itself into clusters or well differentiated pieces. This last statement being a 
consequence of the positive fractal (Hausdorff) dimension associated with a > 0, as 
was mentioned above. 

For a < the situation is qualitatively different. We can distinguish two different 
regimes within non-asymptotically free theories, according as to whether > a > —2 
or —2 > a. In the former case, (r^) diverges as diffusion takes place into larger 
distances, and is only shut-off by the vanishing of a with distance, that is to say, by 
the transition into the classical regime. When a < —2, the coefficient of diffusion 
goes to zero as s increases; this means that new states stationary under diffusion can 
be created and supported at these very deep values of a. When cr > —2 the diffusion 
coefficient diverges for large s and quantum corrections will shut the diffusion process 
off by eventually driving a to zero; what happens is that as s increases, the diffusion 
coefficient at first increases, to then decrease and cut itself off to zero when a goes to 
zero or when the classical regime {h 0) is reached. 

We have seen how quantum corrections lead into diffusion, and in some cases, clus- 
tered states of the quantum field system. To get more information about clustering 
and other structure formation properties, it is convenient to look at the correlation 
function and integral; from an analysis of these objects one can get a better feeling 
of the structures0 supported by the quantum field theory. The correlation integral, 
C{s), gives information on how many correlated pairs there are whose separation 
is less than s. This integral is computed from the correlation function C,{r) by the 
following formula0'0 




V{s) 
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where the region of integration extends over the volume V{s) of radius s. The corre- 
lation function ^(r) is related to S{k), the Fourier transform of the density p{f), via 
the following Fourier transform 



1 + 



where 



S{k) = J , 



The quantity is called the power spectrum, and measures the mean number of 

correlated neighbors in excess of random over a distance of order 1/k; because of this, 
the quantity 1 + ^{f) (once appropriately normalized) gives information on whether 
there is "clustering" ( > 1), "voids" (< 1) or a perfectly random distribution {— 1). 
This is a straightforward program that can be carried out for p(r). One obtains the 
following expressions for 1 + ^(r) and C(r) 



-3+0- 



27r 

C(r) = — • D{a) ■ r 

(7 



2a 



where the coefficient D{a) is 



and 



TT 



D{a) = {AnAy -Att- \B{a)\-T{2-2a) cos-(2a- 1 



TV 



B{a) = r((T-l)cos-(2-a). 

D{a) is shown in Figure 3, and in Figure 4 we show D{(7)/(7, which are important 
in determining the behaviors of 1 + ^(r) and the correlation integral. 



11 



Prom Figure 3 we see that for positive u and between zero and 3/2, the quantum 
field system leads to clustering. For o between 3/2 and 5/2 there is a significant 
qualitative change into a scenario of deep anticlustering ("voids") with a tendency 
for the system to behave near homogeneity when cr 2. When a increases past 
5/2, we enter a new clustering domain. The rest can be "read off' in the same way 
from the figure. The plot of Figure 4 simply confirms the one for D{a), but for the 
coefficient of the correlation integral. We also point out that there is a "quantization" 
of structures, a kind of periodicity, that takes place as a changes and goes through 
some critical values. 

We end by summarizing our results. Quantum field theory is a many body system 
par excellence. Intrinsic to it is the randomness associated with quantum fiuctuations. 
The application of a statistical interpretation to one of the most primitive concepts 
of field theory, that of the induced charge, has interesting and apparently deep con- 
sequences which stem from its scaling properties. In particular, fractal behavior in a 
form related to the Pareto-Levy form of the charge density, leads to systems where 
the maximum entropy is associated with the classical limit of the quantum system. 
In other words, the classical regime is the most stable, and the one preferred by the 
system. 

Like in other many body systems, diffusion is a familiar mechanism for the ap- 
proach to equilibrium (although perhaps not the only one). In concordance to what 
follows from studying the entropy, the virtual cloud diffuses with a coefficient of 
diffusion proportional to h, and thus diffusion of the cloud stops in the classical limit. 

From an analysis of the correlation properties of the virtual cloud system, one 
learns that the system supports the formation of structures where there is a sharp 
transition between phases of anticlustering (void formation), clustering and some 
interspersed quasi-homogeneity. These depend quite clearly on the number of degrees 
of freedom excited from the vacuum into the system; in other words, they depend on 
the distance at which the system is probed. 

In very general terms, we see that irreversibility is contained in quantum field 
systems, and that this irreversibility is encountered as the system size grows. The 
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quantum field system tends to "relax" into a classical system or into clustered struc- 
tures, depending on the scale of the system and the nature of the interactions. The 
ideas developed here may also find application in the study of intermittency, tur- 
bulence, phase transitions, multiparticle physics and other complex phenomena in 
quantum field theory and the early universe. 
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